arXiv:1504.03472v3 [gr-qc] 7 Jul 2016 


Local equilibrium solutions in simple anisotropic cosmological 
models, as described by relativistic fluid dynamics 

Dmitry Shogii|3 and Per Amund Amundsei|3 
Faculty of Science and Technology, 

University of Stavanger, N-fOSd Stavanger, Norway 

Abstract 

We test the physical relevance of the full and the truncated versions of the Israel-Stewart the¬ 
ory of irreversible thermodynamics in a cosmological setting. Using a dynamical systems method, 
we determine the asymptotic future of plane symmetric Blanch! type I spacetimes with a vis¬ 
cous mathematical fluid, keeping track of the magnitude of the relative dissipative fluxes, which 
determines the applicability of the Israel-Stewart theory. We consider the situations where the dis¬ 
sipative mechanisms of shear and bulk viscosity are involved separately and simultaneously. It is 
demonstrated that the only case in the given model when the fluid asymptotically approaches local 
thermal equilibrium, and the underlying assumptions of the Israel-Stewart theory are therefore not 
violated, is that of a dissipative fluid with vanishing bulk viscosity. The truncated Israel-Stewart 
equations for shear viscosity are found to produce solutions which manifest pathological dynamical 
features and, in addition, to be strongly sensitive to the choice of initial conditions. Since these 
features are observed already in the case of an oversimplified mathematical fluid model, we have 
no reason to assume that the truncation of the Israel-Stewart transport equations will produce 
relevant results for physically more realistic fluids. The possible role of bulk and shear viscosity in 
cosmological evolution is also discussed. 
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I. INTRODUCTION 


One of direct consequences of the kinetic theory is that real substances manifest transport 
properties, which leads to existence of dissipative mechanisms like viscosity and thermal 
conductivity. Cosmological fluids should not be an exception. Although perfect fluids in 
local thermal equilibrium (LTE) have been successfully used to model the matter content 
of the universe, it is dissipation that drives the fluid towards equilibrium: a fluid without 
appropriate transport properties will never reach this state. Consequently, dissipation must 
have played an important role in the early universe. 

A brief review of the traditional theories of dissipative thermodynamics can be found 


m 


111 ; 


see references therein for applications to cosmology. At the present moment, the most 


33 , 


often 


advanced of these approaches is the second-order Israel-Stewart (IS) theory 
referred to as transient, or causal, thermodynamics. An alternative formalism proposed 
by Carter j4|, see 3 ^ review, has been shown to be essentially equivalent to the IS 

theory 6|. 

The simplified, or ’’truncated”, version of the IS theory is obtained by omitting certain 
divergence terms in the transport equations of the full theory in a way that does not violate 
the crucial features of causality and stability. As pointed out by Maartens l| , the truncated 
IS theory can be considered as a relativistic generalization of the Navier-Stokes equations 
for a fluid, in the sense that the values of the transport coefficients are equal to their local 
equilibrium values. Due to this property, the truncated IS theory is sometimes treated as 
an independent approach rather than an approximation to the full theory. 


Maartens 3 shown that the truncated IS theory can yield solutions which deviate 
signihcantly from those of the full theory. Furthermore, Zimdahl jsj has demonstated for 
the case of scalar dissipation that the full and truncated versions are equivalent only in 
a few special cases. Recently, Shogin et ah 3 have discovered that the truncation of 
the IS transport equations can produce solutions with pathological dynamical features in 
anisotropic spacetimes. However, despite these results, the truncated theory is still widely 
used in cosmology, see e.g. references in the review by Brevik and Grpn in 101. 

The standard Friedmann-Robertson-Walker cosmological models have been studied using 


the truncated 


11 


■ Il2[ | and the full 3ll2l versions of the IS theory. In particular, the processes 


of inflation and reheating, where dissipative effects play a signihcant role, have been discussed 
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respectively in 



However, the geometry of the spatially homogeneous and isotropic Friedmann-Robertson- 
Walker model does not allow for any dissipative mechanisms different from bulk viscosity. 
Therefore, anisotropic backgrounds must be considered when modelling more realistic dis¬ 
sipation. For an early application of causal thermodynamics to anisotropic cosmological 
models, see Belinskii et ah 1^. Van den Hoogen and Coley 15| have considered Bianchi 
type V cosmological models using the truncated IS theory. Heuristical methods have been 
used to find out that, in contrast to the first-order Eckart theory, the shear viscous stresses 
play a major role in determining the dynamics of the model; namely, depending on the coef- 
hcients of shear viscosity, the spacetime may, or may not, isotropize in the future. Recently, 
Shogin et al. {q] studied Bianchi type IV and V models in the framework both of the full and 
the truncated versions of the IS theory. The existence of effectively anisotropic solutions, as 
well as the importance of shear viscosity, was confirmed; in addition, it was demonstrated 
that even the full IS theory can break down in cosmological applications, as the viscous 
stresses can drive the fluid far away from LTE. 


In this paper, we investigate spatially homogeneous plane symmetric Bianchi type I cos¬ 
mological models. By this choice we intentionally exclude the effects of spatial curvature 
to concentrate on determining the limits of applicability of the IS theory, and, in addition, 
Ending the source of the singular behaviour of the solutions of the system in the case when 
the transport equations are truncated. 


We apply a dynamical systems approach 


16| to describe the future attractors of the 


system of differential equations governing the dynamics of the cosmological model and per¬ 
form a stability analysis of the hxed points of this system. Special attention is paid to the 
late-time behaviour of the relative dissipative fluxes, the magnitude of which represents the 
deviation of the fluid from LTE, and, by this, determines the applicablity of the IS theory 
in this particular case. 


The paper is organized as follows. In section [ITl we briefly discuss the chosen fluid model. 
Then, the Hubble-normalized field equations are written down in section IIIIl The future 
asymptotic states of the cosmological model, as well as the dynamics of the relative dis¬ 
sipative fluxes, are discussed in sections [TVTjVll Conclusions, together with some possible 
directions for the further work are presented in section IVHI 
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II. THE FLUID MODEL 


We consider the cosmological model to incorporate a non-tilted [17|] dissipative fluid, 
neglecting the effects of thermal conductivity. The energy-momentum tensor of the fluid is 
then given by 

Tafi = {p + p + T^)UaUii + {p + + Ta/?, ( 1 ) 

where tt denotes the bulk viscous pressure, Ua is the four-velocity of the fluid, and stands 
for the shear viscous stress tensor, with Ta/su^ = riap] = = 0. We are working in the 

so-called Eckart frame, so the variables p and p refer to the local equilibrium values of energy 
density and pressure of the fluid, and is defined as the particle four-velocity of the fluid. 

Since our purpose is to perform a qualitative investigation, we consider a mathematical 
fluid, for which p and p are connected by a linear barotropic equation of state, widely known 
as the 7 -law: 

p = (7 - l)p, 0 < 7 < 2 . ( 2 ) 


Note that the equations of state for physical radiative fluids are somewhat more compli¬ 


cated, aee e. g. ^ i-oMag temperature and partrcle aurrrber derrsrty at irrdeperrdeat 
variable^. A simple 7 -fluid with 1 < 7 < 4/3 has been used as a very naive model of a 
matter/radiation mixture. 

The dissipative properties of the fluid are described by the IS theory, the transport 
equations of the full versiopj (in a spatially homogeneous background) being formulated 


TqTT -F tt = 


'^2'l~ab T '^ab 


-3CH - -Toll 


C T 


-2paab - ^T2Tah 


To 


T2 V T 


( 3 ) 

( 4 ) 


Here T is the local equilibrium value for the temperature, H is the Hubble rate, and aab 
is the geometric rate of shear tensor. In causal thermodynamics, the relaxation times tq 
and T 2 of the dissipative processes are finite and related to the bulk and the shear viscosity 
coefficients C and p by 


^ A realistic fluid model may, however, lead to simplifications in other aspects of the analysis. 

^ The transport equation describing heat conduction is omitted. 

^ Dots represent derivation with respect to proper time t of observers comoving with the fluid. 
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T2 = 2 ? 7 /? 2 , 


'^0 — C / 3 o 5 


( 5 ) 


where l^o, (^2 > 0 are the thermodynamic coefficients for scalar and tensor contributions to 
the entropy density |l|. 

The transport equations of the truncated IS theory are obtained by dropping the terms 
in the square brackets on the right-hand sides of ([3]) and (jl]), which results in: 


TqTT + 71 = -3CH, 


'^2^ab T '^ab ‘^V^aby 


( 6 ) 

( 7 ) 


while the relations ([5]) still hold. 

To complete the model, the transient and the transport coefficients must be specified. 


Following 


Iq and 9|, we assume a barotropic form of the transport coefficients and the 


relaxation times, the values of the exponents being uniquely determined from dimensional 
analysis: 


C oc 

rj oc 


l3o 

1 

J2 


oc p; 


oc p. 


( 8 ) 


Note that for a realistic fluid these coefficients are more complicated functions of the 


fluid variables and are derived from radiative thermodynamics 


20 


21| . The barotropic 


approximation may be reasonable in some cases, but here we use it purely for the sake of 
technical convenience. 

An artefact of these widely adopted models is the need for a separate equation of state 
for the temperature of the mathematical fluid. For a physical radiative fluid, no special 
’’temperature model” is needed, as temperature is already incorporated in the equations of 


state 


20 


2l| . The two most common choices are 


Q: 


1 . The barotropic temperature model, described by 


T oc p' 


(7-l)/7 


( 9 ) 


and used as an approximation for radiation-dominated mixtures (7 close to 4/3); 
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2. The ideal-gas temperature model, with p = nT, where n is the (conserved) number 
density of the fluid obeying 


VJnu^) = 0 . 


( 10 ) 


In Bianchi type I spacetimes, tensor equation (iTOjl can be rewritten in terms of scalars 
as 

n + 3Hn = 0. (11) 

Note that the physical relevance of this approximation is not proven, with a possible 
exception for ultra-relativistic and highly relativistic ideal gases. 

We use both models in the full version of the IS theory. In the truncated version, how¬ 
ever, a particular equation of state for the temperature will not have any influence on the 
cosmological dynamics, unless a heat-conductive fluid is involved. 

The relative dissipative fluxes caused by the bulk and the shear viscous stresses, respec¬ 
tively, are introduced by: 


Xn = 


X-r = 


{XabT' 


ab\^n 


P 


( 12 ) 


The underlying assumption of the IS theory is that the fluid is close to LTE, which implies 
that these dissipative fluxes are small; 

Xr « 1. (13) 

When condition flT^ is violated, the IS theory breaks down, and the resulting solutions 
cannot be considered physically relevant, although they are mathematically consistent and 
dynamically stable. We do not conjecture that the IS theory can be ’’extrapolated” to 
produce relevant results when applied to fluids far away from LTE; hence, it is extremely 
important to keep track of the dissipative fluxes in all the solutions obtained using the IS 
theory. 

Finally, it is assumed that the mean interaction time ti of the fluid particles is much 
shorter than the characteristic timescale for macroscopic processes, i.e. << for the 
hydrodynamic description to be applicable to the considered matter model. 


III. THE SYSTEM OF EQUATIONS 


We work in the orthornormal frame 


22| and introduce the scale-independent variables 


using the common notations of the held, see e. g. 


16 


23 


24| . The normalized energy 
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density f 2 and the dimensionless geometric shear Tiab are defined by: 


^ = Q, ^ = ^ab = diag(-2S+, S+, S+). (14) 

The bnlk and the shear viscous stresses, respectively, are normalized by: 

^^ = n, ^ = rab = dmg{-2r+,T+,T+). (15) 

For convenience, we consider Bianchi type I spacetimes with planar symmetry, for which 
both Tiab and Tab take the diagonal form. We shall omit the index ” +” in the present paper. 
Then, the Einstein field equations arc: 

S' = (g-2)S + r, (16) 

n' = {2q +2- 37)12 - SB - 2Sr, (17) 


where 

9=-;^-l = -§^-l = 2E"+Q7-l)si + |n (18) 

is the deceleration parameter. In addition, the Hamiltonian constraint binds S and 12 
algebraically: 

1 = S2 + 12. (19) 

The system is completed by the dimensionless transport equations, where the bulk and 

the shear viscosity parameters oq, 5o and 02 , 62 , respectively, are introduced by: 

1 1 

1 1 

1P'J2 

and treated as positive constants. The bulk viscosity parameter oq is in addition restricted 
from above by oq < 87(2 — 7 ), providing that bulk viscous perturbations propagate at finite 
speeds [l|. 

The full transport equations depend on the particular choice of a temperature model. We 
consider the following options: 


0 q 12 , 

02 I 2 , 


1 1 

h'Wo 

1 1 

H 


2t]/3. 


= 6 o\/i 2 , 

- = 62 \/i 2 


( 20 ) 


1. For the barotropic temperature model, the full transport equations are 


H' 


r 


3 l + q 

2 7 

3 l + q 

o - 

2 7 


Bq+T} + 
62 \/H + 


27 - 112' 
27 12 

27 - 112' 
27 12 


n — aol 2 . 


T — 202 ^. 


Primes denote derivation with respect to dimensionless time t, introduced by dt/dt = H. 


( 21 ) 

( 22 ) 
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2. For the ideal-gas temperature model, the full transport equations take the form 


n' = 

r = 


r- n' 

- 6 „x/S) + - 

-(, 2 \/n + k 


n — oohi, 

r-2a2S. 


(23) 

(24) 


3. Finally, the transport equations of the truncated IS theory are: 

n' = 2(g -|- 1) — boy/fi 11 — Oohl, 

T' = 2 (q + 1) — b 2 y/n T — 2a2^- 


(25) 

(26) 


Note that the evolution equations for the temperature T and the particle number density n 
decouple from the main system. Hence, the dynamics of these fluid variables is trivial and 
follows directly from the equation of state and/or particle number conservation. This would 
not be the case in a model incorporating a realistic radiative fluid. 

Note also that the energy density H and the shear stress E are bounded variables, as 
follows from the restriction H > 0 and the Hamiltonian constraint ([19]). However, the state 
space as a whole is unbounded, since no mathematical restriction is imposed on H and T. 
The dimension of the physical state space is three for the fluids with both bulk and shear 
viscosity, and two in case when one of these dissipative mechanisms is not taken into account. 
The state vector then belongs to a subspace of x or x M, respectively. 

In each case, we investigate the full system of equations both analytically and numerically. 
We determine the fixed points of the system and perform an analysis of their local stability in 
the futurtf. This is done by standard analytical methods, see e. g. [^, 3, 261, while multiple 


numerical runs at different sets of model parameters and initial conditions are used to make 
a conjecture about the global attractor of the system. In all the numerical simulations, the 
Hamiltonian constraint is chosen to be initially satished. 

The dissipative fluxes, which determine the applicability of the IS theory, are studied by 
keeping track of the quantities 


X = 


= |7 - l\Xn, 


y = 


r 

n 


= 3|7 — ll^T 


(27) 


® All the fixed points of the considered dynamical systems are not listed in this paper; only those which 
can be locally stable in the future are discussed. 
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IV. SOLUTIONS FOR BULK/SHEAR VISCOUS FLUIDS 


In all the three models under consideration, the only locally stable future attractor is 
described by 


[E, r, Ji, n, (/] = |o,o, i,n,(/|, 


( 28 ) 


with n < 0. This state describes an izotropizing cosmological model, which is dominated 
by the bulk viscous fluid at late times. The sign of q is determined by a model-specihc 
relation between 7 and the bulk viscosity parameters. Thus, depending on this sign, the 
spacetime can end up in a state of decelerated, uniform, or accelerated expansion; the latter 
case represents bulk viscous inflation. Note that anisotropic solutions, which are present in 
anisotropic cosmological models with non-zero spatial curvature 
Bianchi type I backgrounds. 


3 , are not obtained in 


While the numerical runs reveal that all the solutions obtained with the full IS equations 
tend asymptotically to the state fl28|) . the situation is different when the transport equations 
are truncated. For a wide range of values of the viscosity parameters, the solutions behave 
unphysically in the truncated IS theory, running into singularities. We consider this in 
details in section IIV Cl 


For the non-singular solutions, the relative dissipative flux caused by the shear viscous 
stresses decays, while that caused by the bulk viscous stresses tends to a negative constant 
value. In general, this value is not small. Thus, the bulk viscosity can, and generally does, 
prevent the fluid from approaching LTE at late times, which is quite typical for similar 
systems where this dissipative mechanism is involved 9|. 


A. The full transport equations, barotropic temperature 


The constants in 0281) are given by 



760 - \/+ 6700 , 

37 - 2 -F 760 - 7^&o + 67^0 


( 29 ) 
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Calculating the eigenvalues corresponding to the state fl28|) yields 


Ai = 


^2,3 — 


+ 6700 
7 

3n(7 + 1) + 37(7 - 2) - 2762 ± d[3n(7 - 1) + 37(7 - 2) + 2762]^ - 327 


(30) 


^02 


the rational part of A 2,3 being negative. As Ai is negative, the local stability of the future 
attractor requires Re(A 2 , 3 ) < 0. This yields 


02 > ——(sn — 2762)(n + 7 — 2), 

87 

which is fulhlled, since the right-hand side is negative. 


(31) 


B. The full transport equations, ideal-gas temperature 


In this case, the future attractor (1281) is specihed by 

Oq 


n = - 


6 n’ 


2 oo 
2 3 Vo 

and the corresponding eigenvalues are given by 
Ai = —bo, 

3(n + 7 - 2) - 262 ± \/ [3(n + 7 - 2) + 262] ^ - 3202 

A2,3 = -^-• 

Again, Ai is negative; the future stability requires Re(A 2 , 3 ) < 0, which leads to 

3 

0-2 > ^&2(n -F 7 - 2). 


(32) 


(33) 


(34) 


The right-hand side is negative, while 02 > 0. The future attractor is locally stable for the 
whole range of values of the model parameters. 


C. The truncated transport equations 


Now the constants in fl28l) are 

1 r 


n=6 


-37 + bo- a/(37 - boY + 12oc 



37 + &o - V (37 - bo)^ + 12 oo , 


(35) 
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the corresponding eigenvalues being 


Ai = -\/(37 - bo)^ + 12 ao, 


^ 2,3 — 


3 ( 3 n + 37- 2) - 262 ± v'[ 3 (n + 7 + 2 - 262)]^ - 32 o: 


( 36 ) 


The local stability of the future asymptotic state requires Re(A 2 , 3 ) < 0; the solution of this 
inequality yields 


h > ^(37 - 2) + 9n, 

02 > ^ [&2(n + 7 - 2) + 6(n + 7) - 3 (n + 7)^]. 


( 37 ) 


The shear viscosity coefficients turn out to play a crucial role in the truncated IS theory. 
For a wide range of values of these parameters, the requirements (157)) are not fulhlled. In 
this case, no hxed point of the system is locally stable in the future; numerical simulations 
show that all the solutions obtained end up in a singularity. The nature of this singularity 
is that the energy density variable crosses the vacuum boundary (f2 = 0) and becomes 
negative. The transport equations are non-dehned in this region, and the solutions break 
down. This is not found to happen in the full IS theory, but singularities of the same type 
have been reported in more advanced Bianchi models with a dissipative mathematical fluid 
described by the truncated IS theory j^. 

Moreover, the system with the truncated transport equations is found to be sensitive to 
alterations of the initial conditions. That is, it is possible to obtain singular solutions even 
if the viscosity parameters do satisfy the inequalities fl37)) . An example demonstrating this 
feature is shown in Figure [H 

For the cosmological model considered, such singular solutions do not exist in the full IS 
theory. Thus, it is the truncation of the IS equations which allows the energy density to 
manifest unphysical behaviour. 


V. SOLUTIONS FOR BULK VISCOUS FLUIDS WITH VANISHING SHEAR VIS¬ 
COSITY 

In the absence of shear viscosity the future attractor retains the form fl28l) . assumed T = 0. 
The expressions for II and q, obtained in section m also hold in this case. However, the 
instabilities of the truncated IS theory described in section IIV Cl are removed. 
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FIG. 1. The dynamics of ^(t) (left) and S(r) (right) in a Bianchi type I cosmological model with 
a bulk/shear viscous fluid described by the truncated Israel-Stewart theory. The plots illustrate 
that the dynamical behaviour of the solutions depends strongly on the initial conditions. 

At S(0) = — 0 . 2 , n(0) = 0.96 (solid lines), the solution asymptotically approaches the state (|28]l . 
as might be expected. At S(0) = —0.5, fl(0) = 0.75 (dashed lines), the solution ends up in a 
singularity, as the energy density variable crosses the vacuum boundary (Q = 0 ). 

The model parameters are 7 = 6/5, oq = 0.1, 60 = 5.0, 02 = 0.5, 62 = 3.0, so the inequalities (f37l) 
are satisfied. 


Although the future attractor is essentially the same, the dynamics of the approach to it 
is changed, as the corresponding eigenvalues are different from those calculated in section HVl 
It turns out that the shear viscous stresses play only a minor role in the full IS theory, 
provided the bulk viscosity is non-zero; namely, they affect neither the future asymptotic 
state nor its stability. On the other hand, the role of shear viscosity in the truncated version 
of the IS theory is crucial, which is consistent with the results of 15 |. 

The relative dissipative flux caused by the bulk viscous stresses tends to a finite constant 
at late times: in general, the fluid does not asymptotically approach LTE. 


A. The full transport equations, barotropic temperature 


The future asymptotic values of 11 and q are given by (EH]); the corresponding eigenvalues 


are 


Ai = 


y/-k 6700 
7 


A2 = Tn + (7-2)]. 


(38) 
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Since 11 < 0 and 7 < 2 , both eigenvalues are negative reals, and the asymptotic state is 
locally stable in the future. 

B. The full transport equations, ideal-gas temperature 

The values of If and q are provided by fl32l) . the eigenvalues being 

Ai = - 6 o, A 2 = ^[n + (7-2)]. (39) 

Again, the eigenvalues are negative real numbers, and the future attractor is locally stable. 

C. The truncated transport eqnations 

The constants in fl28|) are given by fl35|) . The corresponding eigenvalues are 

Ai = -V( 37 -&o )2 + 12ao, A 2 = ^ [fl + (7 - 2)] . (40) 

In contrast to the case of section IIV Cl the eigenvalues are now real and negative. The 
instability of the future attractor is removed. The singular solutions, which appear in the 
case of nonvanishing shear viscosity, cannot be obtained in this case. 

Hence, the pathological dynamical features of solutions in the truncated IS theory origi¬ 
nate from the evolution of the shear viscous stresses. 

VI. SOLUTIONS FOR SHEAR VISCOUS FLUIDS WITH VANISHING BULK 
VISCOSITY 

The only future stable stationary point of the system is given by 

[S,r,0,g] = [0,0, l,g], (41) 

with 

«= 1(37- 2). (42) 

The asymptotic value of the deceleration parameter depends on 7 only and is positive for 
typical 7 -fluids. Hence, the shear viscosity does not contribute to accelerating the expansion 
of the universe; however, spatial anisotropy is eliminated at late times. 
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For the full IS transport equations, the eigenvalues are either negative reals or complex 
conjugates with a negative real part; this provides the local stability of the future attractor. 
For the truncated transport equations, the situation is different: the negativity of Re(Ai^ 2 ) 
can be violated, which makes the state fHT]) unstable in the future. The corresponding 
solutions are found to end up in a singularity, which is discussed in section IVI Bl below. 

For all the non-singular solutions, the relative dissipative fluxes caused by the shear 
viscous stresses decay in the future, and the fluid approaches LTE at late times. This is the 
only case in the considered cosmological model when the underlying assumption (IT^ of the 
IS theory is not violated and the solutions of the system can be considered fully reasonable. 


A. The full transport equations 

The eigenvalues corresponding to the state fl4ip are the same for both temperature models: 

^ _ (37 - 262 - 6) ± V(37 + 262 - 6)2 - 32a2 

^1,2 - -^ 

where the rational part is negative. The local stability requirement Re(Ai^ 2 ) < 0 yields 

a2>^&2(7-2), (44) 

which is fulhlled automatically by 02 > 0 and 7 < 2 . 


B. The truncated transport equations 


For the truncated transport equations, the eigenvalues are given by 

_ (97 - 262 - 6 ) ± V (37 - 262 + 6)2 - 3202 
^1.2 --^-• 


(45) 


An algebraic analysis of the local stability requirement Re(Ai^ 2 ) < 0 results in the following 
system of inequalities: 


62 >5(37-2). 

“2 > j (7 - 2)(62 -37). 


(46) 


If the values of 02 and 62 do not satisfy these conditions, the corresponding solutions behave 
unphysically and end up in a singularity of the same nature as described in section IIV Cl 
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Moreover, similarly to the case with nonvanishing bulk and shear viscosity, the stability 
of the future attractor depends crucially on the initial conditions; it is possible to obtain 
singular solutions, even if the inequalities are satished. As might be expected, this is 
not a property of the system with the full transport equations. 

This conhrms our assumption that it is the evolution equations for the shear viscous 
stresses that lead to unacceptable properties of the solutions obtained using the truncated 
IS theory. 


VII. CONCLUSIONS 


We have used the dynamical systems approach to investigate the future attractors and 
their stability conditions for viscous mathematical fluids in Bianchi type I spacetimes. We 
have studied the properties of cosmological solutions obtained with the full and the truncated 
versions of the IS theory, having used two simple temperature models in the full version. 
Also, we have determined the asymptotic future of the relative dissipative fluxes to hnd out 
when the near-equilibrium conditions are violated and the IS theory breaks down. 

All the solutions obtained using the full IS theory are found to be non-singular in future 
and stable under alterations of the initial conditions. On the contrary, the truncated IS 
transport equations, if applied to a fluid with nonvanishing shear viscosity, allow the energy 
density to cross the vacuum boundary, which leads to unphysical, singular behaviour of the 
solutions already in the simplest anisotropic spacetimes. 

Another pathological feature of the truncated IS equations discovered in the present 
work is the extreme sensitivity of the resulting system to the choice of the initial conditions. 
When the truncated IS theory is applied to a 7 -fluid with nonvanishing shear viscosity, the 
solutions can run into a singularity even if there exists a stable future asymptotic state. 

The solutions of the full IS theory describe an izotropizing universe dominated by the 
dissipative fluid in the asymptotic future. The shear viscous stresses, if present, decay at late 
times, while the bulk viscous stress, if present, freezes into a negative constant value. In these 
solutions, the bulk viscosity eliminates spatial anisotropy and, in addition, can accelerate 
the expansion of the universe. The shear viscosity does not contribute to accelerating the 
expansion, but effectively eliminates the Bianchi type I anisotropy even when bulk viscosity 


is zero. Anisotropic solutions, which exist in more complicated spacetimes 


15|, cannot 
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be obtained in Bianchi type I cosmological models. 


The full IS theory provides a completely reasonable description only for shear viscous 
mathematical fluids with vanishing bulk viscosity: the relative dissipative flux caused by 
the shear viscous stresses decays exponentially in the future, and the fluid is driven towards 
LTE. 


The full IS theory yields stable solutions for bulk viscous fluids. However, these solutions 
are not fully consistent with the underlying assumptions of the IS theory, since generally the 
fluid neither is close to LTE during its evolution nor approaches thermal equilibrium in the 
asymptotic future. The near-equilibrium conditions can be violated by the bulk viscosity: 
the corresponding relative dissipative flux is asymptotically a negative constant, which is 
not small in general. This departure from LTE is, of course, a physical result, which may 
depend on the fluid model adopted, and per se not a failure of the IS theory. 


Large deviations from LTE caused by the shear viscous stresses, which have been discov¬ 
ered in more general cosmological models {q], are not present in Bianchi type I spacetimes. 


The future asymptotic states of the Bianchi type I spacetimes are of similar form for the 
two temperature models considered in the full IS theory. Still, the solutions can behave 
differently in the future; for example, these models under the same initial conditions can 
result in opposite signs of the asymptotic value of the deceleration parameter q. Also, 
the dynamical character of the solutions (monotoneous or oscillatory behaviour) can differ 
between them. 


The results obtained suggest that in solving real physical problems, the full IS theory 
should be preferred over its truncated version. For a mathematical fluid with barotropic 
transport coefficients, the bulk viscosity creates hnite, non-decaying relative dissipative 
fluxes, which lead to a possible breakdown of the full IS theory; hence, consistent non¬ 
linear thermodynamical theories, which can describe fluids substantially away from LTE, are 
potentially important. However, for realistic radiative fluids with transient and transport 
coefficients derived from the kinetic theory, the dynamics of the variables may be essentially 
different. We leave considering physical fluids in anisotropic cosmological backgrounds for 
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the further investigations. 
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